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Abstract 

In this paper we study a family of operators dependent on a small 
i parameter e > 0, which arise in a problem in fluid mechanics. We 

show that the spectra of these operators converge to N as e — > 0, even 
though, for fixed e > 0, the eigenvalue asymptotics are quadratic. 
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1 Introduction 

> 

■ In a recent paper [1] Benilov, O'Brien and Sazonov argued that the equation 

oo ■ 

q: i-Hf (i; 

^ ■ approximates the evolution of a liquid film inside a rotating horizontal cylin 

O ■ der, where H is the closure of the operator Hq on L^(— vr, vr) defined by 

W)W^4(sm(.)|)+| (2) 

for any sufficiently small fixed e > and all / G Dom (i/o) = Cper([~^' ^])- 
They also made several conjectures, based on non-rigorous asymptotic and 
numerical analysis, including that the spectrum of H is purely imaginary and 
consists of eigenvalues which accumulate at ±ioo, and that the eigenvalues 
converge to iZ as e — 0. Weir proved these conjectures, except for the 
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convergence of the eigenvalues, in We prove the remaining conjecture 

in this paper. 

Davies showed in [1] that, for < e > 2, every A G C is an eigenvalue of 
H . For e < 2 he showed that —iH has compact resolvent by considering the 
unitarily equivalent operator A on /^(Z) defined by 

{Av)n = ^n{n - l)vn-i - ^n{n + l)vn+i + nVn (3) 

for all V e Dom(A) = {v e : Av G /^(Z)}. Here A = J^-\-iH)J^, 

where JF : 7r,7r) — > Z^(Z) is the Fourier transform. If JF/ = v then 

{vn)n£Z are the Fourier coefficients of /. This result was achieved by obtaining 
sharp bounds on the rate of decay of eigenvectors and resolvent kernels, and 
by determining the precise domains of the operators involved. He also showed 
that 

A = A_ © © A+, (4) 

where A_ and A^ are the restrictions of A to /^(Z_) and P{Z^) respectively, 
and that A- is unitarily equivalent to —A+. Since the resolvent is compact 
and the adjoint has the same eigenvalues, the spectrum of —iH consists 
entirely of eigenvalues. 

Weir proved in PJ that these eigenvalues, if they exist, must all be real. 
Boulton, Levitin and Marietta subsequently proved in a recent paper [S] 
that a wider class of operators possess only real eigenvalues. However, they 
did not prove that any non-zero eigenvalues exist for these operators, nor 
that their spectra are real. In [3], Weir proved rigorously that —iH has 
infinitely many eigenvalues, all of multiplicity one, which accumulate at ±oo 
by showing that the eigenvalues of ^4+ correspond to those of a self-adjoint 
operator with compact resolvent. This operator on L'^{{0,l),2ws{x)dx) 
is defined as the closure of the operator given by 

iLef){x) = -'-W,{xr\pj'nx) (5) 

for all f e{f e C°°([0, 1]) : /(O) = 0}, where 

We{x) = X-\l - xY/%1 + xy'/' , 

Pe{x) = (l-x)^+^/"(l + a;)i-^/^ 

It was argued in [6] that the distribution of the eigenvalues, if they exist, 
should be quadratic, but no rigorous bounds were given. By analysing the 
self-adjoint operator. Weir proved rigorously in [3] that A„ ~ £7r^n^/9~^ as 
n oo for a certain explicit constant (3. 
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The conjecture in P that A^^^ — > as e — >• is supported by numerical 
evidence in the same paper, and also in [HI SI E]. In this paper we prove 
the conjecture rigorously. In Section [2] we apply unitary transformations 
to obtain a family of operators on the same space. These operators are 
invertible with Hilbert-Schmidt inverses (Theorem \2.2\\ . We then identify a 
differential operator with spectrum N (Theorem 13.31 Corollary 13.41) . which 
is unitarily equivalent to an operator whose inverse has an integral kernel 
which is the pointwise limit of the integral kernels of the aforementioned 
family of inverses (Theorem 13.61) . Finally we show that we actually have 
norm convergence (Theorem 14. 2p and use the variational method to show 
that this implies convergence of the eigenvalues, i.e. X^^n n as e —>■ for 
all n G N; see Theorem 14. 3[ 



2 Unitary transformations 

Since we are interested in the limit as e ^ 0, we assume for the rest of the pa- 
per that < e < 1. In order to obtain convergence of operators in some sense, 
we need a family of operators on the same space. To this end we apply uni- 
tary transformations to to obtain a family of operators on L'^{{0, 1/e), ds) 
and then extend each operator to an operator on L^{{0, oo), ds). 

Lemma 2.1 The operator Lg, defined by is unitarily equivalent to an 
operator on L'^{{0, l/e),Ws{s)ds) such that 

iL,g){s) = -w,{s)-\p,g'ns) (6) 

where 

p,{s) = {l-esy+^/'il + esY'^/' 
for all s G (0, 1/5). Moreover, is invertible with inverse Rs given by 

iRef)is)= G,{s,t)f{t)We{t)dt (7) 

for all f G L'^{{0, 1/e), We(s)ds), where 



7e(s) = / ps{u) Mm 
Jo 

1 ( n+es\^/' 



2 VI -es 



(9) 
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and 



Ge{s,t) 



%{s) zfO<s<t<l/e 
jelt) ifO<t<s<l/e. 



(10) 



Proof We define an operator U : -^^^((0, 1), 2ws{x)dx) -^^^((0, l/e),We{s)ds) 
by 

iUf){s) = f{es) (11) 

for all / G L'^{{0,l),2ws{x)dx), s G (0, 1/e). It is easy to show that U is 
a unitary operator. We define = ULg^U^^. A simple calculation shows 
that ^ holds. The analogous result about the inverse of was proven in 
[3], so the result about follows from the unitary equivalence of the two 
operators. It only remains for us to calculate 



le[.S) 



e Jo 



e Jo 1 



- 

1 



e JO (1 — u) 



1~J 

1 + u 



du 

l/e-l 



du 



u 



u 



1/e 



l + es \ 
l-esj 



J 

l/e 



Theorem 2.2 The operator is unitarily equivalent to an operator on 
L^((0, l/e), ds) such that 



for all f G L'^{{0, l/e), ds), where 



l/e 



K,{s,t)f{t)dt 



l/2e 



K.[s,t) 



my" 

l/e 



l-st \ 
1+etJ 



l/2e 
l/2e 



(12) 



ifO<s<t 
ifO<t<s. 



(13) 
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Moreover, if we define on L^((0, oo), ds) by 

{NJ){s)= / K,{s,t)f{t)dt (14) 
Jo 

for all f G -Z^^((0, oo), ds), where 

I U otherwise, 

then iVg has the same non-zero eigenvalues as , and each non-zero eigen- 
value has the same multiplicity with respect to each operator. 

Proof We define a unitary operator Je : L^((0, 1/e), ds) ^ L^((0, 1/e), We(s)ds) 
by 

{J,f){s)=w,{s)-'/^f{s) (16) 
and then put = J'^L^Je. We then have = J^^R^Je, so 

{M;'f){s)= f''K,{s,t)f{t)dt (17) 
Jo 

for all / e Iv2((0, l/e), ds), where 

K,{s,t) = Ms)'^'Gs{s,t)ws{t)'/' 

j w,(s)V2^,(s)w},(t)i/2 if < s < t < l/e, 

\ We{s)^/^-fe{t)We{&^^ if < t < S < l/e. 

Substituting in expression for 7^, we obtain ( fT3l) . 

There is no unitary equivalence between and N^, but it is easy to 
see that they have the same eigenvalues, and that the non-zero eigenvalues 
have the same multiplicities with respect to each operator. Indeed, 

N, = © 0, 

where is the zero operator acting on L^((l/£:, 00), ds). The theorem now 
follows from the unitary equivalence established in Lemma 12.11 | 

3 The limit operator 

In this section we consider the operator Lq on H = L^((0, 00), u'o(s)ds), 
where 

wois) = Mmweis) = 2s-^e-^\ (18) 
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defined on some suitable domain, which we identify below, by 



(Lo/)(.) = -Wo{s)-\pof'ns), (19) 

where po{s) = lim^^o p^{s) = e~^'^. We shall show that Lq has a self-adjoint 
extension, which is invertible and whose inverse is Hilbert-Schmidt. We then 
identify the spectrum of Lq and the integral kernel of its inverse. 

Definition Let V denote the set of all polynomials on (0, oo) and sV denote 
those elements of V which have constant term zero. 

Definition Let Dom (Lq) be the set of twice differentiable functions g : 
(0, oo) C such that lims^o+ 9{s) = 0, lim sup^^Q_,_ \9'{s) \ < oo, lims^oo Q~'^g{s) 
0, lims^oo e~'^(7'(s) = and Lq/ G TC. 

Remark It is easy to show that Dom (Lq) C Ti. We may therefore define 
Lq on this domain by equation ( fT9l) . Note also that sV C Dom {Lq). 

Lemma 3.1 The setV is dense in L^((0, oo), 2se~^'^ds). 

Proof Let g E L'^{(0, oo), 2se~^''ds) be such that g -LV and define 

for all s e R. We need to prove that g = 0, or equivalently that / = 0. We 
first prove that / G L'^[TV) and the Fourier transform / of / has an analytic 
continuation to the set {z : l^^l < 1}. We then prove that / is zero on an 
interval containing zero, and hence is identically zero. The invertibility of 
the Fourier transform then implies that / = 0. 
Firstly 



/ \f{s)fds< sup 2se~^' r \g{s)f2se-^'ds < 
"'R. se(o,oo) "'0 



oo 



so / G L^(R). If 6 < 1 then a similar estimate shows that e^''''/ G L^(R). 
Theorem IX. 13 of [7] implies that / has an analytic continuation to {z : 
\^z\ < 1}. 

In order to evaluate / in a neighbourhood of zero, we first approximate 
g-i^x polynomials: 



e 



^ n' 

n=0 "'• 



< E 



n=N 



6 



= \x\ 



< 



I l^lel^ 



E 

n=0 



N\ - ni 



for all a; e R, ^ e R, e N. We now choose /3 G (0, 1). If 5 > is sufficiently 
small then e^^+'^^l^l/ e L'^{R). Also e"-^!^! e ^^(R), so e^l^l/ e Li(R). If 
1^1 < P/2 then 



/ e-^^^/(x)d, 



< 



R 



<-L 



x 

N-l 

E 

n=0 



ni 



f{x)dx 



+ 



R 



n=0 



ni 



m 



\f{x)\dx + 



N-l 

E 



-i^x) 



< sup 

a;6(0,oo) 



I l^lcl^ 



AT! 



n=0 

/ |/(a;)|e^l"lda; 

^R 



-g{x)2xe dx 



for all A^ e N. The final integral is finite since e^^^^f E L^(R). We put 
h^,N{x) = -l^LMLed^l-^)!^! for all x e (0, 1) and note that this is a smooth, 
positive function of x. As x — > or x ^ oo, 0, so the supremum of h 

must be obtained at a local maximum in the interval (0, oo). The only zero 
of h' is at N/{P - |^|) so 

sup h{x) = hiN/iP-m 



N 



N 







as AT ^ oo for all ^ e [-/3/2, /3/2]. Hence /(O = for all ^ e [-/3/2, /3/2]. 



Corollary 3.2 The set sV is dense in H. 
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Proof It is easy to show that the operator U : L^((0, cxd), 2se '^ds) — > 7i 
defined by 

{Uf){s)^sf{s) 
for all s e (0, oo) is unitary and that UV — sV. | 

Theorem 3.3 The operator Lq is symmetric and there is a complete or- 
thonormal sequence {e„}„gN in H such that, for each n e e„ e Dom (Lq) 
and LoCn — ncn- 



Proof For all f,ge Dom (Lq) 

{Lof,g) = - / {pofns)g{s)ds 
Jo 



po(0)/'(0)^(0) - Imi po{n)f{n)g{n) + / po{s)f{s)g'{s)ds 



= \\mpo{n)f{n)g'{n)-po{0)f{0)g'{0)- I f{s){pog'y{s)ds 
= if: Log) 
since /(O) = ^(0) = 0, 



lim Po{n)f'{n)g{n) = ( lim e ^f'{n) I ( lim e ^g{n) I = 



and similarly \imn^(X}Po{n)f{n)g'{n) = 0. Therefore Lq is symmetric. 
For all n e N, define a„^j. recursively for r = 1, . . . , n by 

an,n = 1 (20) 
vir -\- 1) 

= "7^7 7«n,r+i for r = 1, . . . , n - 1 (21) 

and put 

fn{s)^Y.an,rS' (22) 

for all s e (0, cxd), n G N. Then, for each n e fn ^ Dom (Lq) and 

= -|r (^^) + ^/'(^) 

r=2 ^ r=l 

r(r + lKr+i\ 

= j'^Qn.r ^ J * ""^"'"^ 

= w/n(s) 
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for all s G (0, oo). Since the fn are eigenvectors corresponding to dis- 
tinct eigenvalues and Lq is symmetric, they are orthogonal. Putting e„ = 
fn ||/riir\ we obtain an orthonormal sequence such that Loe„ = ne„. Each 
e„ is a polynomial of degree n, so lin{e„}„gN = sV. It now follows from 
Corollary 13.21 that {e„}„gN is complete. | 

The situation above can be seen as an exceptional case of that for the asso- 
ciated Laguerre polynomials {Ll^^}n=o (see 0). If a > — 1 these polynomials 
are a complete orthogonal set of eigenvectors, corresponding to eigenvalues 
0,1,2,... respectively, of the differential equation 

xf" + (a + 1 - x)f' + A/ = 

on L'^{{0, oo),x°'e~^dx). After a change of variables, the eigenvalue equation 
for Lq becomes the exceptional case a = —1. The singularity at the origin 
requires special treatment, but otherwise the treatment of this case is the 
same as for a > — 1. 

Corollary 3.4 The operator Lq is essentially self-adjoint, and the spectrum 
of its closure is precisely N. 

Proof The result follows immediately from the theorem and Lemma 1.2.2 
of [g. I 

Theorem 3.5 The operator Lq is invertible and its inverse R is a Hilbert- 
Schmidt operator. Moreover, 

oo 



{Rof){s)= / GQ{s,t)f{t)wQ{t)dt (23) 
Jo 

for all / G 7Y, where 

r (, f\-l ^ois) zfO<s<t 
°^ ' ^" 1 7o(t) ^fO<t<s 

and 

7o(s) = / pQ{u)~^du 
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Proof Let Rq be defined by equation fl23l) . Then Rq is Hilbert-Schmidt since 



oo / roo 



\J0 



\Go{s,t)f wo{t)dt] wo{s)ds 



oo / rco 



oo / /"OO 



Go{s,t)f wo{t)dt) wo{s)ds 

2 



< 2 



< 



oo / poo 



e'^' - 1 



e^^-l 



e-2*dt 

2 



t-^e-2*dt ) s-^e-^Ms 



l„-2s 



s-\-''ds 



Ads+ / 



Let / G sP. Then / G i^^((0, oo), iL'o(s)ds). Hence it is easy to show that 
Rf is absolutely continuous and 



{RfYis) 



i{s)f{t)womt 



Po[s) 



fit)woit)dt. 



We now show that Rf G Dom(Lo). Since /, Wq and p^^ are smooth on 
(0,oo), (Rf) is twice differentiable (indeed, it is also smooth). We calculate 



limsup|(i?/)(s)| <limsup7o(s) / \fit)\woit)dt = (25) 



and 



limsup|(i?/)'(s)| < limsuppo(s)~W \f{t)\wo{t)dt 

coo 

\f{t)\Wo{t)dt<^. 

Since / is a polynomial, there exists a constant c > such that |/(s)| < ce'^ 
for all s G (0, oo). Let 5 > be given. Then, for all s > c/6, 



< e-' 

< e~' 



ho{t)fit)woit)\dt + e 



hois)fit)woit)\dt 



/ e2*|/(t)|rV2*dt + e- / e'^\fit)\ 

Jo Js 



t-'e-^'dt 



c/S 



\fit)\t-'dt + e 



6e'dt + s'^e' I ce-'dt 



c/5 



c/5 



\f{t)\r^dt + 6 + cs-^ 6 
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as s — oo. Since 5 > is arbitrary, e (s) — > as s ^ oo. Also 

1"°° . , 2 



|/(t)|-e-2*dt 

< 2s- vy 



roo 



= 2cs-^ 

as s — s> oo. Therefore Rf G Dom (Lq), provided L^Rf G 7^, and 
{L,Rf){s) = {LoRf){s) 



oo 



= -Ms)''-^J^ f{t)wo{t)dt 
= fis). 

Since / G 7i, we indeed have L^Rf G 7i. 

We have proven that L^Rf = f for all / G sP. Let / G 7i. By Corollary 
13.21 there is a sequence (/„) in such that /„ — > / as n ^ oo. By the 
above, we have Rfn — > Rf and L^Rfn = /„ / as ^ oo. Since Lq 
is closed, this implies that Rf G Dom (^Lq) and LoRf = f. Finally, Lq is 

injective since ^ Spec (^Loj . Hence L^RLof = L^f imphes RLof = f for 

all fen. I 

Theorem 3.6 The operator Lq is unitarily equivalent to an operator Mq on 
L^((0, oo), ds). Moreover, 



(Mo-V)(s)= / Ko{s,t)f{t)dt (26) 
Jo 

for all f G L'^{{0, oo), ds), where 

Ko{s,t) = wo{sy^^Go{s,t)wo{ty/^ 

j (St)-l/2g-s |g2s _ 1} e"* %fO<s<t 

j (st)-i/2g-^{e2* - l}e-* z/0<t<s. 
Proof We define a unitary operator 

Jo : L2((0, oo), ds) L\{0, oo), Wo(s)ds) 

by 

{Jof){s)=w,{s)-'/'f{s) (27) 
and put Mq = Jq^LqJq. Then Mq^ = Jq^RqJq and fpUl) is immediate. | 
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4 Convergence of the eigenvalues 

Clearly K£{s, t) — > Ko{s, t) as e ^ for each s, t e (0, oo). In this section we 
shall show that Kgr — ^ Kq in L^-norm as e — > 0, and hence that A^^ — > Mq^ 
in Hilbert-Schmidt norm as £ — > 0. We use this fact to prove that the 
eigenvalues of converge to those of Lq as £ — > 0. 

Lemma 4.1 IfO<s<t< 1/e then 



s-t 



Proof We first note that 



(28) 



.og((af)"^'(sf)"^' 



2s 



{I0g(l+£S) -I0g(l-£S)} 



and hence 



+ -{log(l-£t)-log(l + £t)} 



oo ^2k 



< S-t 



I /2s 



1-es) \l + st) 



(29) 



for < s < t < l/£. For such s, t, 



1— £6- 

1+es 



1/2e 



{^ f -1 

1—es 



1+st 



< 

< e 



l+es 
1—es 



l/2e 



\ 1+et J 



l/2e 



{e2*-l}e-* + e-*-*. 



Theorem 4.2 linie^o 



HS 



0. 



Proof Using the symmetry of and Kq, it is sufficient to show that 



Js 



Ks-K, 



dtds 



(30) 
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as e — * 0. If h\og2 < s < t < l/e then, by Lemma F4.lt 



St 



Ke-Ko 



l—es 
1+es 



l/2e 



1+£S 

l—es 



l/e 



1-et 
1+st 



l/2e 



|g2s _ 



< max<ie ^*,e ^'^ 



< e-^'^je^^-ll'e-^*. 
If t > l/e then i^^j = 0, so the above bound still holds. Since 

poo pco , ^2 

/ / (st)-V2^|e2^-ire"2*dMs 



< 



ilog2 



{e--l] 



e-2*dMs 



2 Jilos2 



< 

- 2 7iiog2 

< CXD, 



we may use the Lebesgue Dominated Convergence Theorem to prove that 

2 



oo roc 



ilog2 Js 



dtds 



(31) 



as £ ^ 0. 

If < s < 1 and < e < 1 then 



log 



l — es 



2E 



k=0 



2k + 1 



< 2E 



oo g2k+l 



k=0 



log 



2A; + 1 



1 + s 



so 



Also 



l+es\^/' ^ , 

-. < log 

1 — esj 

1 -£X\l/2e 



1 - s 

1 + s 



1 - s 



1 + ex 



< e" 



(32) 
(33) 
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if < a; < Hence, for < s < i log 2, s < t < 

2 



< St 



< max 



< max |e-2- (i±f - l)' e'^', e'^' {e^^ - 1}' e-2*| 



As before, this bound holds for all t such that t > s, since — if t > 1/e. 
Since we have 



5-V^^{e2^-l}'ds 



< 



2 Jo 

< oo 



2 /"OO 



e-2*dMs 



1 /•5l°g2 



and 



1 + s 



-1] e-^'dtds 



< c/' / e-2^e-2*dMs 

C /■l'*^ -4,, 
= 2/. 



< 00 



for some constant c, we may use the Lebesgue Dominated Convergence The- 
orem to prove that 



ilog2 roc 



Jo Js 



dtds 



(34) 



as £ — > 0. 



We now use standard variational methods to deduce the convergence of 
the eigenvalues from the norm convergence of the resolvents. 

Theorem 4.3 For each e G (0, 1), n e let X^^n be the nth eigenvalue of 
L^. Then Xs^n — > n as £ — >• 0. 
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Proof If S* is a finite-dimensional subspace of L'^{{0, oo), ds) then we define 



fieiS) = sup{((/-iV,)/,/):/G Sand 11/11=1}, (35) 

/io(S) = sup{((/-5o-^)/,/):/G5and 11/11 = l}, (36) 

fis^n = inf{/i£(5') : S C L^((0, oo), ds) and dim(S') = n}, (37) 

/io,n = inf{/io('S') : S C L^((0, oo), ds) and dim(S') = n}. (38) 

For each e G (0, 1), A^e is self-adjoint since it is a Hilbert-Schmidt operator 
with a symmetric integral kernel. Also Spec (L^) C [1, oo) and hence it 
follows from Theorem 12.21 that J — A^^ is a non-negative self-adjoint operator 
with essential spectrum {1} and eigenvalues {1 — '■ n G N} C [0, 1). It 

follows from Corollary 13.41 and Theorem 13.61 that / — Mj"^ is a non-negative 
self-adjoint operator with essential spectrum {1} and eigenvalues {1 — 1/n : 
n G N} C [0,1). 

By Theorem 4.5.2 of [9], fie,n = 1 — 1/Ae,n and /io,n = l — l/n. Let 5 > be 
given. Then Theorem 14.21 implies that there exists r/ > such that, whenever 
< e <r], 

((/ - Mo"^)/, /) - 5 < ((/ - N,)f, /) < ((/ - Afo"')/, /) + 6 (39) 
for all / G -^^^((0, oo), ds) such that ||/|| = 1. For such e, this implies that 

f^eiS) - 5 < floiS) < fleiS) + 5 (40) 

for all finite-dimensional subspaces S of L'^{{0, oo), ds) and hence that 

Ate,n -S < yUo,n < fJ'e,n + S (41) 

for all n G N. Therefore, for all n G N, fi^^n — ^ /^o,n as e — > and hence 
Ae,n ^ n as £ ^ 0. I 
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